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Ultrafast  optical  control  of  entanglement  between 
two  quantum-dot  spins 

Danny  Kim,  Samuel  G.  Carter,  Alex  Greilich,  Allan  S.  Bracker  and  Daniel  Gammon* 

The  interaction  between  two  quantum  bits  enables  the  creation  of  entanglement,  the  two-particle  correlations  that  are  at  the 
heart  of  quantum  information  science.  In  semiconductor  quantum  dots,  much  work  has  focused  on  demonstrating  control  over 
single  spin  qubits  using  optical  techniques.  However,  optical  control  of  two  spin  qubits  remains  a  major  challenge  for  scaling 
to  a  fully  fledged  quantum-information  platform.  Here,  we  combine  advances  in  vertically  stacked  quantum  dots  with  ultrafast 
laser  techniques  to  achieve  optical  control  of  the  entangled  state  of  two  electron  spins.  Each  electron  is  in  a  separate  InAs 
quantum  dot,  and  the  spins  interact  through  tunnelling,  where  the  tunnelling  rate  determines  how  rapidly  entangling  operations 
can  be  carried  out.  We  achieve  two-qubit  gates  with  an  interaction  rate  of  30  GHz,  more  than  an  order  of  magnitude  faster 
than  demonstrated  in  any  other  system  so  far.  These  results  demonstrate  the  viability  and  advantages  of  optically  controlled 
quantum-dot  spins  for  multi-qubit  systems. 


Semiconductor  quantum  dots  (QDs)  were  among  the  first 
candidates  proposed  for  solid-state  qubits1.  Self-assembled 
InAs  QDs  are  a  particularly  versatile  platform,  because  they  are 
epitaxially  grown  in  a  semiconductor  wafer  and  can  be  fabricated 
into  a  monolithic  architecture  containing  both  electronic2  and 
photonic3  circuit  elements.  Individual  QDs  themselves  can  be 
organized  into  more  complex  ‘molecules’  in  one,  two  and  three 
dimensions4.  With  these  engineering  advantages,  one  can  anticipate 
an  entire  quantum  network5  with  the  scalability  and  stability  of  a 
solid-state  system6. 

The  elementary  optical  excitation  of  a  QD,  the  exciton,  has  a 
resonance  frequency  in  the  optical  regime,  giving  QDs  a  great  speed 
advantage  over  nuclear  spins  (radiofrequency)  or  electron  spins 
(microwaves) .  With  the  giant  optical  dipole  of  a  semiconductor  QD, 
quantum  operations  can  be  carried  out  at  a  terahertz  rate  or  faster7. 
Coherent  manipulations  of  exciton  qubits  were  some  of  the  earliest 
quantum  gate  demonstrations  in  the  solid  state8-11 . 

Unfortunately,  the  exciton  lifetime  is  less  than  a  nanosecond 
in  a  QD,  so  in  practice  they  are  inadequate  as  qubits.  This  prob¬ 
lem  is  solved  by  charging  the  QD  with  a  single  electron,  and 
using  the  long-lived  electron  spin  as  the  qubit.  The  electron  spin 
lives  for  milliseconds12  or  longer  and  has  a  coherence  time  of 
microseconds13-16.  The  exciton  now  acts  as  an  auxiliary  state  for  ul¬ 
trafast  conversion  between  optical  coherence  and  spin  coherence17. 

The  opportunity  created  by  this  optically  controlled  spin 
paradigm  has  produced  a  burst  of  recent  activity  in  single-qubit 
operations  with  QDs,  including  single-spin  initialization18-20,  non¬ 
destructive  readout20,21  and  fast  spin  manipulation22-26.  However, 
the  power  of  quantum  information  originates  from  entanglement 
of  multiple  qubits,  and  so  far  all  of  these  efforts  are  limited  to  a  single 
isolated  qubit.  Creating  and  controlling  entanglement  requires  a 
well-controlled  spin  interaction  between  QDs,  and  quantum  optical 
techniques  that  can  address  individual  QDs  within  a  pair.  The 
interaction  needs  to  be  strong  enough  to  maintain  the  speed 
advantage  of  optically  controlled  QDs,  but  not  so  strong  that 
individual  control  of  dots  is  impractical. 

Here,  we  go  beyond  one  QD  qubit,  demonstrating  the 
fabrication  and  optical  control  of  a  two-qubit  system  consisting 
of  two  electron  spins,  each  on  separate  tunnel-coupled  QDs. 


We  establish  concurrently  the  three  important  prerequisites 
for  two-qubit  entanglement  control:  (1)  initialization  of  the 
two-electron  spin  state  through  optical  pumping;  (2)  single-qubit 
gates  using  short  laser  pulses;  (3)  and  two-qubit  gates  using  longer 
laser  pulses.  To  achieve  single -qubit  gates  in  a  multiple -qubit 
system,  the  key  is  to  use  a  pulse  so  short  that  the  two  dots  do  not 
have  time  to  interact  while  a  single  spin  is  optically  rotated.  For 
two-qubit  gates  we  use  pulses  that  are  longer  than  the  interaction 
time  between  the  dots.  Combining  all  of  these  techniques,  we  obtain 
optical  control  of  two  QD  qubits  on  a  timescale  that  is  much  faster 
than  any  other  candidate  for  quantum  computing27-29. 

Tunnel-coupled  QDs 

Controlled  entanglement  of  two  spins  requires  an  interaction 
between  the  spins,  the  strength  of  which  determines  how  fast  the 
two-qubit  gate  can  be  carried  out.  A  number  of  entanglement 
mechanisms  have  been  proposed  for  optically  driven  QD  spins, 
including  coupling  through  photons  and  optical  cavity  modes17,30,31 , 
and  through  Coulomb  and/or  tunnelling  interactions  between 
neighbouring  dots32-36.  Here  we  employ  the  kinetic  exchange 
interaction  based  on  coherent  tunnelling.  Tunnelling  provides  the 
largest  interaction  rate,  and  because  of  recent  advances  in  the 
growth  and  spectroscopy  of  vertically  stacked  QDs  (refs  37-40) 
can  be  precisely  controlled  through  a  combination  of  quantum 
size  engineering  and  applied  electric  field  bias.  This  interaction 
has  also  been  used  in  electrically  defined  QDs27,41,  although  the 
interaction  is  much  stronger  in  our  self- assembled  QD  system — 
30  GHz  as  compared  with  ~2GHz  (ref.  27),  allowing  an  order  of 
magnitude  faster  gate  speeds. 

For  this  study,  two  vertically  stacked  self- assembled  InAs  QDs 
were  grown  with  a  thin  tunnel  barrier  of  GaAs/AlGaAs  such  that  the 
electrons  can  coherently  tunnel  between  the  dots  (Fig.  1).  We  grew 
the  two  QDs  with  different  thicknesses  so  that  they  have  different 
optical  transition  energies.  As  a  result  each  is  optically  addressable 
with  a  resonant  laser  frequency.  The  QDs  are  incorporated  into  a 
Schottky  diode  so  that  by  adjusting  the  voltage  bias  each  QD  is 
charged  with  a  single  electron. 

The  calculated  Hund-Mulliken  model  of  the  2e  energy 
level  system42  is  shown  in  Fig.  2.  Using  optical  transmission 
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Figure  1 1  Sample  structure.  Schematic  of  the  Schottky  diode  structure  by 
which  two  electrons  are  charged  into  separate  QDs.  The  spins  in  the  two 
QDs  are  coupled  by  coherent  tunnelling  through  a  9  nm 
GaAs/Alo.sGaojAs/GaAs  (3/3/3)  barrier.  Laser  fields  are  used  to 
initialize,  measure  or  rotate  the  spins  through  a  real  or  virtual  electron-hole 
pair  (exciton)  resonantly  excited  in  the  top  dot. 


spectroscopy,  the  2e  charge  state  was  found  to  be  stable  for  a 
bias  range  of  100  mV  centred  at  about  0.46  V.  Beyond  this  range 
1  and  3  electrons  are  more  stable  in  the  QD  pair.  Within  this 
2e  stability  range  the  lowest  energy  configuration  consists  of  one 
electron  residing  in  each  dot.  As  a  result  of  the  Pauli  exclusion 
principle,  the  two  electrons  can  tunnel  between  QDs  only  if  they 
have  an  antisymmetric  spin  state.  The  spin  states  are  split  into 
degenerate  triplets  (T0  =  ltl)  +  llt)>  T-  =  I>I4)>  T+  =  |tt))  and 
a  lower-energy  singlet  (S  =  lt^)  —  l^t))  with  a  kinetic  exchange 
splitting  Z\ee  that  is  a  result  of  tunnelling  between  the  dots43. 
The  magnitude  of  the  exchange  energy  determines  the  interaction 
strength  between  dots  and  ultimately  how  fast  the  entanglement  can 
be  controlled.  We  found  that  an  exchange  energy  of  30  GHz  enables 
both  ultrafast  one-  and  two-qubit  gates,  as  we  will  show. 

The  spin  state  is  optically  controlled  through  an  exciton  state 
(X2-),  in  which  an  electron-hole  pair  is  created  in  the  top  dot. 
An  exciton  in  the  smaller  bottom  dot  is  tens  of  millielectronvolts 
higher  in  energy.  Tunnelling  between  the  two  dots  is  rendered 
negligible  in  the  optically  excited  state  at  the  same  bias  (0.46  V) 
owing  to  additional  Coulomb  interaction  energies42,44.  We  start 
at  zero  magnetic  field,  in  which  case  the  exciton  state  is  fourfold 
degenerate,  with  the  states  labelled  by  their  total  spin  projection  mz 
in  Fig.  2.  The  optical  selection  rules  at  zero  field  (blue  lines)  show 
that  both  S  and  T0  couple  to  the  same  optically  excited  state,  forming 
a  A -level  diagram,  as  required  for  optical  spin  control. 

Two-qubit  initialization:  Optical  spin  pumping 

Optical  spin  initialization  of  the  two-electron  (2c)  system  is 
demonstrated  by  optical  pumping.  The  two  spectral  lines  in  Fig.  3a 
arise  from  excitation  of  the  2e  singlet  and  triplet  states  to  the 
common  exciton.  The  splitting  between  the  lines  gives  a  direct 
measure  of  the  kinetic  exchange  energy  (Z\ee  ~  125  peV  =  30  GHz). 
Figure  3  a  shows  bleaching  of  the  singlet  and  triplet  absorption  lines 
in  the  bias  range  between  0.45  and  0.48  V,  that  is  the  centre  of 
the  2e  charge  stability  plateau.  This  is  a  well-known  signature  of 
optical  spin  pumping  in  a  A -level  diagram18-20,  where  the  spin  is 
excited  out  of  one  spin  state  and  becomes  shelved  in  the  other 
after  a  few  recombination  cycles.  The  pumping  is  efficient  (at  least 
95%  fidelity)  as  long  as  the  optical  cycle  is  faster  than  relaxation 


0.3  0.4  0.5  0.6  0.7  0.8 

Applied  bias  (V) 


Figure  2  |  Energy  levels  of  the  2e  system.  Hund-Mulliken  model  of  the  2e 
system  showing  the  energy  levels  of  the  ground  and  lowest  optically 
excited  states,  and  their  respective  tunnelling  anticrossings.  The  red  arrow 
shows  approximately  the  bias  in  which  these  experiments  were  conducted. 
Inset:  Corresponding  energy  level  diagram  for  spin  substates  (not  to  scale). 
The  X2-  excited  trion  states  are  labelled  by  their  total  spin  projection  mz. 
The  selection  rules  at  zero  magnetic  field  are  shown  by  blue  arrows 
(circularly  polarized).  In  a  transverse  magnetic  field  additional  transitions, 
shown  by  red  arrows,  become  fully  allowed,  and  the  blue  and  red 
transitions  are  linearly  polarized  and  orthogonal. 

between  the  S  and  T  ground  states.  Near  the  edges  of  the  stability 
plateau  (white  dashed  lines  in  Fig.  3a),  spin-flip  cotunnelling  to 
the  doping  layer  results  in  efficient  ground-state  relaxation,  and 
no  net  pumping  is  possible.  The  spin  pumping  can  be  defeated  by 
repumping  on  the  opposite  branch  of  the  A -system  using  a  second 
laser  as  shown  in  the  two  linescans  in  Fig.  3c. 

The  selection  rules  of  Fig.  2  (inset)  indicate  that  two  of  the  triplet 
sublevels  (T+  and  T_)  are  not  part  of  a  A -diagram  and  therefore 
should  not  pump  at  zero  magnetic  field.  However,  as  shown  by 
the  experiment,  the  triplet  absorption  signal  is  almost  completely 
suppressed  through  optical  pumping.  This  implies  that  there  is 
weak  mixing  of  the  triplets  through  the  hyperfine  interaction 
and/or  heavy-hole/light-hole  mixing  in  the  exciton  states.  These 
interactions  turn  on  additional  weak  optical  transitions  and  allow 
full  initialization  into  the  singlet  state. 

Although  spin  initialization  through  optical  pumping  may  seem 
routine,  in  this  case  incoherent  pumping  is  in  fact  generating  a 
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Figure  3  |  Optical  spin  pumping  of  the  two-electron  state,  a,  Intensity  plot 
of  the  single-laser  transmission  spectrum  as  a  function  of  bias.  The  zero  in 
the  energy  scale  is  at  1.3  eV.  The  dashed  horizontal  lines  indicate  the  edges 
of  the  two-electron  stability  plateau.  Inset:  A-diagram  with  the  singlet  and 
triplet  coupling  to  the  same  X2-  state,  b,  Single-laser  linescan  spectra  at 
two  bias  values:  one  in  the  cotunnelling  regime  that  shows  triplet  and 
singlet  peaks  (0.51  V)  and  one  in  the  pumping  regime  where  the  peaks  are 
gone  (0.46  V).  c,  These  traces  show  that  the  transmission  signals  at  0.46  V 
are  restored  by  tuning  the  frequency  of  a  second  repumping  laser  to  the 
opposite  arm  of  the  A-system. 

particular  entangled  spin  state,  S  =  1 1|)  —  |||).  This  entanglement 
is  possible  owing  to  the  exchange  interaction.  It  may  seem  strange  at 
first  to  consider  the  individual  spin  states  as  qubits  and  the  singlet  an 
entangled  state  because  the  exchange  interaction  is  always  coupling 
the  two  spins.  Nevertheless,  we  can  independently  manipulate  the 
individual  spin  states  using  optical  pulses  that  act  faster  than  this 
interaction.  This  allows  us  to  transform  the  singlet  state  into  a 
variety  of  other  entangled  states. 

Entanglement  control  using  pairs  of  single-qubit  gates 

For  single-qubit  gates  we  use  short,  circularly  polarized  pulses 
(Dps)  with  a  bandwidth  (146 peV)  larger  than  the  exchange 
splitting  but  much  less  than  the  splitting  between  excitons  in  the 
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top  and  bottom  dots  (Fig.  4a).  As  the  pulses  are  shorter  than  the 
exchange  interaction  time  and  because  the  exciton  is  localized  in 
the  top  dot,  the  pulse  acts  only  on  the  spin  in  this  dot.  The  pump 
laser  pulses  are  detuned  below  the  triplet  transition  by  ~270  peV, 
allowing  a  single-spin  rotation  with  only  virtual  occupation  of  the 
X2-  exciton  state22,24,26.  This  single-qubit  gate  induces  a  relative 
phase  between  the  up  and  down  spin  states  in  just  the  top  dot — a 
rotation  of  this  spin  about  the  optical  axis.  As  shown  previously  for  a 
single  electron  spin,  a  circularly  polarized  Jt-rotation  takes  1 1)  —  1 1) 
into  ||)  —  eZJt|t)  =  ||)  +  ||)>  (refs  22,24-26,45).  Similarly,  in  this 
two-electron  case,  where  the  pulse  acts  only  on  the  top  dot,  the 
jt-pulse  rotates  S=  It  l)  — lit)  into  T0  =  1 1 1)  + 1 1 1)- 

By  applying  two  such  single-spin  rotations  separated  by  a  delay 
time  during  which  the  2e  spin  state  precesses,  we  can  control  the  2e 
state.  Initialization  into  the  singlet  state  again  is  achieved  with  the 
continuous -wave  (cw)  probe  laser  tuned  to  the  triplet  transition, 
which  also  acts  as  the  measurement.  At  zero  magnetic  field,  we 
operate  in  the  S  and  T0  subspace,  which  can  be  visualized  as  a 
Bloch  sphere  as  shown  in  Fig.  4b.  The  first  pulse  rotates  the  2e 
state  from  S  to  a  superposition  with  T0,  with  a  rotation  angle  that 
depends  on  the  laser  power.  For  example,  an  effective  jr/2-pulse 
rotates  the  Bloch  vector  to  the  equator  as  shown  in  Fig.  4b.  Between 
pulses,  the  state  precesses  around  the  Bloch  sphere  at  the  frequency 
of  the  exchange  splitting  (Z\ee).  The  second  pulse  drives  the  state 
either  up  or  down  depending  on  the  phase  of  the  superposition 
state.  The  projection  onto  T0  determines  the  population  and  the 
corresponding  signal  of  the  cw  probe. 

Figure  4c  (upper  curves)  shows  the  oscillations  in  the  signal  as 
a  function  of  pulse  delay,  called  Ramsey  fringes.  The  frequency  of 
these  oscillations  (30  GHz)  corresponds  to  the  exchange  energy  Z\ee, 
which  is  a  function  of  bias  voltage.  The  observed  frequencies  as  a 
function  of  bias  (red  dots  in  Fig.  4d)  are  in  excellent  agreement  with 
the  measured  singlet-triplet  splitting  (red  dashed  curve  in  Fig.  4d). 
The  oscillations  have  a  decay  time  of 400-700  ps,  depending  on  bias, 
although  this  can  be  substantially  increased.  The  main  damping 
source  is  found  to  be  fluctuations  in  electric  field  that  produce 
fluctuations  in  the  exchange  energy.  This  is  seen  in  Fig.  4d,  where 
the  bias  dependence  of  the  decay  has  been  fitted  to  the  derivative 
of  the  exchange  energy  (<5Z\ee/<5V)  multiplied  by  the  voltage 
fluctuation  amplitude  of  3.8  mV,  taken  as  the  fitting  parameter. 
These  fluctuations  probably  arise  from  background  charges  in 
the  sample.  This  contribution  can  be  removed  by  designing  the 
molecule  such  that  <$Z\ee/<5V  =  0  within  the  stability  range  of 
the  2e  configuration.  This  would  be  accomplished  by  shifting 
the  stability  region  of  the  2e  configuration  in  Fig.  2  towards  the 
right  to  ~0.6V.  Additional  dephasing  contributions  due  to  the 
cw  initialization/measurement  laser  that  is  persistently  on  can 
be  eliminated  by  turning  this  laser  off  during  the  precession  in 
between  pulses16.  Hyperfine-induced  coupling  to  the  T+/T_  states 
is  expected  to  lead  to  dephasing  on  a  longer  timescale,  although  a 
number  of  recent  results  on  single  QDs  indicate  that  even  this  effect 
can  be  substantially  reduced13-15. 

The  amplitude  of  the  oscillations  as  a  function  of  pulse  intensity 
is  plotted  in  Fig.  4e.  The  maximum  amplitude  occurs  for  two 
tc/2 -pulses,  but  decreases  to  nearly  zero  for  two  Jt-pulses,  which 
should  drive  the  system  back  to  the  initial  S  state,  regardless  of  pulse 
delay.  These  oscillations  as  a  function  of  laser  power  correspond 
to  spin-Rabi  oscillations,  in  which  the  Bloch  vector  can  be  driven 
all  the  way  around  the  Bloch  sphere.  The  amplitude  dependence 
on  the  pulse  intensity  is  not  linear  as  expected  from  the  adiabatic 
approximation.  Instead  the  dependence  is  sublinear  (P0'7),  which 
indicates  that  the  detuning  is  insufficient  for  adiabatic  elimination 
of  the  exciton  state24.  The  combination  of  laser  pulse  area  and 
delay  time  permits  the  controlled  generation  of  any  entangled 
superposition  state  within  this  part  of  the  2e  spin  phase  space, 
that  is  a|  t|)  +  /3|||). 
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Figure  4  |  Coherent  control  of  the  2e  spin  state  using  two  13  ps  pulses  with  a  variable  time  delay,  a,  Single-qubit  gate  schematic  where  a  short  pulse  is 
used,  with  a  bandwidth  greater  than  the  singlet-triplet  splitting,  b,  Bloch  sphere  representation  of  spin  control  within  the  S-To  subspace.  The  spin  state  is 
initialized  to  the  south  pole  (S  =  Iff)  -  )).  The  first  pulse  coherently  drives  the  Rabi  vector  up,  after  which  it  precesses  around  the  sphere.  The  second 

pulse  drives  the  vector  either  up  or  down  depending  on  the  time  separation  between  pulses  and  thereby  leads  to  oscillations  in  the  measurement,  c,  Upper 
curve:  the  resulting  Ramsey  interference  oscillations  in  the  To  transmission  signal  as  a  function  of  the  pulse  separation  time,  with  zero  magnetic  field. 
Lower  curve:  Ramsey  interference  obtained  at  0.2  T  showing  two  frequency  components,  comprising  the  faster  S-To  splitting  (30  GHz)  and  the  slower 
T+/T_  splitting  (2.2  GHz).  The  latter  frequency  is  at  twice  the  single-electron  Zeeman  splitting  (ge  =  -0.4).  d,  Measured  Ramsey  fringe  frequency  and 
decay  time  as  a  function  of  bias  voltage  at  0  T.  The  dashed  line  for  the  fringe  frequency  is  the  measured  singlet-triplet  splitting  (Z\ee)  from  Fig.  3a.  The 
dashed  line  for  the  decay  time  is  a  calculation  of  the  expected  dephasing  time  due  to  voltage  fluctuations  of  Z\ee  as  discussed  in  the  text,  e,  Fringe 
amplitude  versus  pulse  intensity  showing  Rabi  oscillations  of  the  spin  state.  The  solid  line  is  a  fit  giving  a  single-qubit  rotation  fidelity  of  98%.  The  error 
bars  are  obtained  from  the  fits  to  the  Ramsey  fringes  at  each  power  and  represent  one  standard  error. 


Optical  control  of  the  full  2e  spin  state  is  possible  using  a 
transverse  magnetic  field  that  allows  Raman  transitions  to  the  T+ 
and  T_  states46  (see  Fig.  2,  inset).  Figure  4c  (lower  curves)  shows 
Ramsey  fringes  in  a  small  transverse  magnetic  field  of  0.2  T,  which 
now  show  frequency  components  associated  with  T+  and  T_.  In 
a  transverse  magnetic  field,  the  spin  eigenstates  are  in  a  different 
basis  with  spins  oriented  along  the  magnetic  field  instead  of  the 
optical  axis.  In  the  optical  basis,  a  short  pulse  rotates  S  ->  T0.  It 
can  be  seen  by  rewriting  S  and  T0  in  the  field  basis  that  a  short 
pulse  now  rotates  S  T+  —  T_ .  The  superposition  of  T+  and 
T_,  Itt)  —  |||)>  precesses  at  twice  the  single- electron  Zeeman 
frequency.  Calculated  Ramsey  fringes  in  this  magnetic  field  are 
shown  in  Fig.  4c  and  reproduce  the  experimental  result  very  well.  In 
this  model,  the  populations  of  the  T0  and  T±  states  after  both  pulses 
are  (1/2)W2[1  —  cos(&>+  —  co_)At]  and  W (1  —  W )[1  +  cos&>±At], 
respectively.  The  frequencies  co±  represent  the  S-T±  splittings,  W 
is  related  to  the  single  pulse  area  O  through  W  —  sin2  ((9/2),  and 
At  is  the  delay  between  pulses.  We  have  observed  that  by  using  a 


Ji-pulse,  full  rotation  to  this  T+ /T_  subspace  is  possible  where  then 
only  these  slower  oscillations  at  co+  —  co _  are  observed.  These  results 
show  that  we  have  extended  the  subspace  from  S-T0  to  include 
the  T+  and  T_  states. 

Two-qubit  phase  gate 

It  may  seem  surprising  in  these  Ramsey  fringe  experiments  that  we 
can  control  the  entangled  state  using  only  single-qubit  rotations. 
In  reality  the  exchange  interaction  that  causes  precession  on  the 
Bloch  sphere  acts  as  a  two-qubit  gate,  with  the  timing  between 
pulses  effectively  controlling  the  interaction  time.  A  direct  way 
to  control  the  2e  state  is  with  longer,  narrowband  pulses  that 
act  as  two-qubit  phase  gates.  In  analogy  with  pulse  rotations 
developed  for  single  QDs  (refs  26,45),  a  narrowband  pulse  near¬ 
resonant  with  an  energy  eigenstate  that  drives  the  system  up 
to  an  excited  state  and  then  back  down  (a  2jr-pulse)  induces  a 
phase  change.  For  example,  an  optical  2jt-pulse  driving  the  singlet 
state  gives  |tl)  —  lit)  =>  e^(ltl)  —  llt))>  where  0  depends  on 
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Figure  5  |  Two-qubit  phase  gate,  a,  Two-qubit  gate  schematic  in  which  a  long  pulse  acts  only  on  one  of  the  ground  states,  b,  With  a  pulsed  pump  laser 
(~150  ps  pulse  width)  on  the  singlet,  optical  Rabi  oscillations  are  observed  with  a  period  that  scales  with  the  square  root  of  the  pump  intensity,  c,  Pulse 
sequence  for  a  two-qubit  phase  gate  experiment,  where  the  second  short  pulse  has  variable  time  delay,  d,  Ramsey  interference  oscillations  with  and 
without  a  two-qubit  control  pulse  near-resonant  with  the  triplet.  The  fidelity  of  the  two-qubit  phase  gate  is  about  80%. 


pulse  detuning25,26,45.  When  this  gate  with  0  =  jc/2  (often  called 
square-root-of-swap)  acts  on  a  product  state,  lt^)>  the  result  is  an 
entangled  state,  | f  i)  —  i || t ) • 

To  obtain  this  two-qubit  phase  gate,  we  first  demonstrate  optical 
Rabi  oscillations  (Fig.  5b).  A  pulsed  laser  with  a  relatively  narrow 
bandwidth  of  12  peV  (~150ps  pulse  width),  much  less  than  the 
singlet-triplet  splitting,  is  used  to  coherently  drive  the  singlet 
transition.  The  cw  probe  is  tuned  to  the  triplet  transition.  As 
the  pulse  intensity  increases,  the  triplet  transmission  signal  shows 
Rabi  oscillations  that  are  periodic  with  the  field  amplitude  (square 
root  of  the  average  laser  power).  The  absorption  signal  does  not 
return  to  zero  owing  to  recombination  during  these  long  pulses, 
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as  reproduced  by  a  simulation  using  optical  Bloch  equations  in  a 
three-level  model  (red  curve  in  Fig.  5b). 

The  two -qubit  phase  gate  is  demonstrated  by  placing  a  2  jt -pulse, 
in  this  case  near-resonant  with  the  triplet  transition,  after  the  first 
short  pulse  of  a  Ramsey  fringe  experiment  (see  Fig.  5a, c).  Pulse  1 
rotates  the  Bloch  vector  up  near  the  equator,  where  it  precesses, 
and  the  longer  control  pulse  changes  the  phase.  Pulse  2  is  used 
to  measure  these  dynamics  as  a  function  of  time  by  rotating  the 
vector  up  or  down,  depending  on  its  phase.  In  Fig.  5d,  the  oscillating 
signal  goes  nearly  to  zero  during  the  control  pulse,  as  the  triplet 
population  is  excited  to  the  exciton.  As  the  exciton  is  driven  back 
down  to  the  triplet,  the  oscillations  return  with  a  phase  change  of 
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145°.  We  have  been  able  to  vary  this  phase  between  — 180°  and  180° 
by  changing  the  pulse  detuning. 

Outlook 

We  have  demonstrated  ultrafast  optical  control  of  two  electron 
spins  in  two  separate  QDs  using  optical  initialization,  single-qubit 
gates  with  short  pulses,  and  two-qubit  gates  with  longer  pulses 
or  through  precession  in  the  exchange  field.  Local  entanglement 
between  the  two  spins  is  inferred  from  the  coherent  evolution  of 
superposition  states  as  measured  in  Ramsey  fringes.  The  fringes 
indicate  that  we  have  generated  a  nearly  pure  quantum  state  of  two 
spins,  which  cannot  be  written  as  a  product  state.  A  direct  measure 
of  entanglement  could  be  carried  out  by  probing  correlations  in  the 
individual  spin  states  on  a  timescale  faster  than  the  interaction. 

One  might  ask  whether  spin  qubits  based  on  the  singlet  and 
triplets  or  the  individual  spin  states  are  most  natural.  These 
pictures  are  formally  equivalent,  and  both  are  useful.  Our  main 
motivation  for  treating  the  qubit  system  as  two  individual  spins 
is  that  we  can  manipulate  the  spins  independently  with  ultrafast 
pulses.  However,  we  find  that  the  simplest  explanation  of  the 
Ramsey  fringes  is  provided  in  the  S-T0  basis  using  the  Bloch  sphere 
illustration  of  Fig.  4a,  which  naturally  accounts  for  precession 
in  the  exchange  field.  This  same  duality  is  present  in  many 
systems,  including  NMR-based  quantum  computing  in  molecules47 
and  quantum  transport  in  electrically  defined  QDs  (ref.  27,48). 
Moreover,  in  quantum  transport  studies,  it  has  been  pointed  out 
that  a  S-T0  qubit  is  of  special  interest  because  of  its  weaker 
sensitivity  to  nuclear  fluctuations. 

An  important  remaining  question  is  what  operations  are 
sufficient  and  most  convenient  to  obtain  universal  control,  allowing 
us  to  span  the  entire  phase  space.  A  well-known  route  to  universal 
control  is  single-qubit  gates  for  each  QD  and  an  entangling 
two-qubit  gate49.  Taking  this  approach  would  require  single-qubit 
gates  about  an  arbitrary  axis  for  both  QDs.  Extending  short-pulse 
spin  rotations  to  the  bottom  QD  should  be  straightforward,  and 
rotation  about  an  arbitrary  axis  has  been  demonstrated  in  single 
QDs  (ref.  24).  A  potential  challenge  in  our  system  is  to  carry 
out  single- qubit  rotations  about  the  magnetic  field  axis  without 
a  simultaneous  two -qubit  interaction.  This  is  a  consequence 
of  having  the  exchange  interaction  in  the  ground  state,  which 
continually  couples  the  qubits  as  they  are  also  precessing  about 
the  magnetic  field.  Ultimately  it  may  be  useful  to  adapt  techniques 
from  NMR,  such  as  refocusing  pulses  to  effectively  turn  off 
the  exchange  interaction47.  Decoupling  the  spins  in  the  ground 
state  could  also  be  achieved  through  quantum  state  engineering. 
For  example,  we  may  engineer  the  coupling  of  the  spin  qubits 
to  occur  only  in  the  exciton  states20,40,43  through  variations  in 
structural  design.  This  would  turn  on  the  interaction  only  when 
required,  and  the  ultrafast  optical  techniques  demonstrated  here 
would  still  be  applicable. 

We  emphasize  that  one  of  the  main  advantages  of  this  optically 
functional  QD  system  is  its  natural  coupling  to  photons,  which  can 
act  as  flying  qubits  for  quantum  communication.  For  example,  the 
spin  entanglement  in  the  dots  could  be  transferred  to  entanglement 
of  emitted  photons50.  Perhaps  most  exciting  is  the  potential  for 
further  scaling  using  a  quantum  network,  an  approach  to  quantum 
computing  or  communication  in  which  a  set  of  nodes  are  entangled 
through  quantum  channels5,30.  These  locally  entangled  spin  qubits 
are  ideal  as  a  node  in  a  distributed  quantum  network  connected 
through  photons6.  For  error  suppression  it  is  necessary  to  have 
within  each  node  at  least  two  qubits  in  which  fast  entanglement 
control  is  possible,  as  demonstrated  here6. 

Methods 

Sample  structure.  Two  vertically  stacked  InAs  dots  were  grown  by  molecular 
beam  epitaxy  on  an  n+ -doped  GaAs  substrate  (see  Fig.  1).  The  QDs  are  embedded 


in  a  GaAs  Schottky  diode,  which  allows  for  deterministic  charging  and  for  tuning 
of  energy  levels  through  the  quantum-confined  Stark  effect.  The  sample  consists 
of  an  n-type  Si- doped  buffer  layer,  a  40  nm  GaAs  spacer  barrier  and  the  two 
InAs  dot  layers  separated  by  a  9  nm  GaAs/AlGaAs/GaAs  barrier.  A  280  nm  GaAs 
capping  layer  with  a  25  nm  AlGaAs  current  blocking  layer  was  included  10  nm 
from  the  surface.  The  top  dot  was  grown  thicker  such  that  it  would  have  the 
lower  optical  transition,  and  such  that  electrons  would  tunnel  but  holes  would 
not44.  The  nominal  dot  thicknesses  are  2.6  and  3.2  nm  for  the  bottom  and  top 
dot,  respectively.  To  obtain  an  appropriate  tunnelling  rate  we  increased  the  height 
of  the  barrier  by  incorporating  3  nm  of  Alo.3Gao.7As  in  between  two  3  nm  layers 
of  GaAs;  giving  a  total  barrier  thickness  of  only  9  nm.  The  actual  anticrossing 
energy  is  2 te  =  700  peV  as  measured  from  photoluminescence.  With  pure  GaAs, 
a  20  nm  barrier  would  be  required  to  obtain  the  same  tunnelling  rate.  The 
energy  between  the  singlet  and  triplet  (exchange  energy)  at  the  bias  where  the  2e 
configuration  is  stable  was  ~  125  peV.  The  sample’s  optical  transition  spectrum  and 
charging  behaviour  as  a  function  of  bias  were  determined  by  photoluminescence 
as  discussed  in  earlier  studies44. 

Measurement  method.  For  single-molecule  studies,  1  pm  apertures  were  patterned 
in  a  120-nm-thick  A1  layer  using  electron-beam  lithography.  A  5  nm  layer  of  Ti  was 
deposited  before  the  Al.  The  sample  was  placed  in  a  He-flow  cryostat  at  5  K.  Piezo 
positioners  were  used  for  positioning.  A  0.68-numerical-aperture  aspheric  lens  was 
used  to  focus  the  laser  to  a  diameter  of  ~2  pm  on  the  aperture,  and  also  to  collect 
the  light  for  the  photoluminescence  measurements.  Light  transmitted  through  the 
aperture  was  focused  onto  an  avalanche  photodiode.  The  transitions  were  probed 
by  means  of  Stark-shift  modulation  spectroscopy  using  narrow-linewidth  lasers 
(nanoelectronvolt  range)  and  lock-in  techniques.  The  sample  was  modulated  with 
a  100 -mV- amplitude  square-wave  voltage  at  10  kHz  superimposed  on  the  fixed  d.c. 
voltage.  A  300  ms  time  constant  was  used  for  all  of  the  experiments. 

The  single  cw  laser  experiments  were  carried  out  with  linearly  polarized 
light.  For  the  Rabi  oscillation  and  Ramsey  fringe  experiments,  one  or  two 
titanium:sapphire  lasers  were  used,  with  circular  polarizations  opposite  to  that 
of  the  cw  probe.  The  pump  laser(s)  for  all  of  the  multi-laser  experiments  were 
rejected  using  a  polarization  analyser.  For  the  Ramsey  fringe  experiments,  short 
pulses  ( 13  ps)  were  detuned  below  the  triplet  transition  by  270  jieV,  and  long  pulses 
(~  150  ps),  if  present,  were  tuned  resonant  with  the  triplet. 

Simulations.  The  Rabi  oscillations  shown  in  Fig.  5b  were  simulated  with  a 
three-level  A-system,  in  which  S  and  T0  couple  to  the  same  trion  state  (one  of 
the  mz  =  ±  1  states)  with  equal  dipole  moments.  The  three-level  density  matrix 
equations  are  solved  numerically  in  response  to  a  120  ps  Gaussian  pulse  resonant 
with  S.  A  recombination  time  of  500  ps  is  used  (equally  to  both  S  and  T0),  and 
recombination  during  the  pulse  prevents  the  population  from  being  driven 
fully  back  down  to  S. 

Simulations  of  the  Ramsey  fringe  experiment  were  also  carried  out,  with  12  ps 
pulses  detuned  300  peV  from  the  centre  of  the  transitions.  These  pulses  rotated 
the  Bloch  vector  from  S  towards  T0,  with  little  trion  population  remaining  after 
the  pulse.  When  the  delay  between  pulses  was  varied,  the  T0  population  due  to 
the  pulses  oscillated  as  observed  experimentally.  However,  the  pulses  were  not 
short  enough  to  completely  rotate  the  Bloch  vector  all  the  way  to  the  T0  state.  We 
estimate  that  the  series  of  two  pulses  used  in  these  experiments  were  able  to  give  a 
maximum  of  65%  population  in  the  T0  state.  We  have  demonstrated  a  full  rotation 
to  the  T0  state  by  using  shorter  3  ps  pulses. 

Fidelity  estimates.  The  initialization  fidelity  of  95%  is  estimated  from  the 
bleaching  of  the  triplet  line  in  the  single-laser  differential  transmission  scans  in 
Fig.  3b.  At  0.46  V,  where  pumping  occurs,  any  residual  absorption  is  clearly  within 
the  noise,  so  we  use  the  noise  level  as  an  upper  bound  on  the  residual  signal.  We 
compare  this  noise  level  with  the  triplet  signal  at  0.51  V,  where  pumping  is  inhibited 
by  rapid  cotunnelling.  The  fidelity  is  then  Fvumv  =  1  —  01 /a0  =  0.95,  where  a  is  the 
noise  level  at  0.46  V  and  a0  is  the  signal  level  at  0.5 1 V. 

We  estimate  the  fidelity  of  single-qubit  gates  using  data  from  Fig.  4e,  which 
plots  the  Ramsey  fringe  amplitude  as  a  function  of  pulse  power.  The  oscillations 
show  decay  that  we  attribute  to  pulse-induced  decoherence.  From  the  fit  to  a 
decaying  cosine,  we  determine  the  expected  amplitude  for  two  Jt/2-pulses  with  no 
dephasing,  Afit,  and  compare  this  with  the  maximum  signal  at  0.3  mW,  Ameas.  If 
each  pulse  reduces  the  Bloch  vector  by  a  factor  D,  then  D2  =  Ameas/Afit  =  0.935, 
giving  a  gate  fidelity  of  F  =  (1  +D)/2  =  0.98  ±  0.01  (ref.  24).  Nominally,  this  is  the 
gate  fidelity  for  a  Jt/2-pulse,  but  these  13  ps  pulses  are  not  short  enough  to  rotate 
the  Bloch  vector  all  the  way  to  the  equator.  From  simulations,  pulses  that  should 
give  a  Jt-pulse  for  a  negligible  exchange  splitting  only  rotate  the  Bloch  vector  by  a 
polar  angle  ~ji/3.  This  effect  is  not  included  in  the  fidelity  calculation  as  it  does 
not  degrade  the  coherence  of  the  state. 

The  fidelity  of  the  two -qubit  narrowband  pulse  is  estimated  from  the  ratio 
of  the  Ramsey  fringe  amplitude  before  the  pulse  (33  ps)  and  after  (443  ps).  This 
gives  D  =  0.61  (F  =  0.80),  which,  surprisingly,  is  better  than  that  expected  from 
dephasing  during  this  time  period,  exp(— 410/500)  =  0.44.  We  speculate  that  the 
pulse  is  partially  reversing  the  inhomogeneous  dephasing  from  voltage  fluctuations, 
but  this  effect  is  not  yet  understood. 
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